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Abstract 

In Refs. [1] and [2], calculation of effective resistances (two-point resistances) on 
distance-regular networks was investigated, where in the first paper, the calculation was 
based on the stratification of the network and Stieltjes function associated with the net- 
work, whereas in the latter one a recursive formula for effective resistances was given 
based on the Christoffel-Darboux identity. In this paper, evaluation of effective resis- 
tances on more general networks called pseudo-distance-regular networks [3] or QD type 
networks [4] is investigated, where we use the stratification of these networks and show 
that the effective resistances between a given node such as a and all of the nodes /3 be- 
longing to the same stratum with respect to a {R^^(m) , (3 belonging to the m-th stratum 
with respect to the a) are the same. Then, based on the spectral techniques, an analyti- 
cal formula for effective resistances Ra/3(m) such that = L^^ (those nodes a, P of the 
network such that the network is symmetric with respect to them) is given in terms of the 
first and second orthogonal polynomials associated with the network, where is the 
pseudo-inverse of the Laplacian of the network. Prom the fact that in distance-regular 
networks, = L^^ is satisfied for all nodes a, /? of the network, the effective resistances 

g{m) for m — 1,2, d (d is diameter of the network which is the same as the number 
of strata) are calculated directly, by using the given formula. 

Keywords:effective resistance, pseudo-distance-regular networks, stratifi- 
cation, spectral distribution 

PACs Index: 01.55. -^b, 02.10.Yn 
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1 Introduction 

The study of electric networks was formulated by Kirchhoff [5] more than 150 years ago as 
an instance of a linear analysis. Our starting point is along the same line by considering the 
Laplacian matrix associated with a network. The Laplacian is a matrix whose off-diagonal en- 
tries are the conductances connecting pairs of nodes. Just as in graph theory where everything 
about a graph is described by its adjacency matrix (whose elements is 1 if two vertices are 
connected and otherwise), everything about an electric network is described by its Lapla- 
cian. Recently, the authors have given a method for calculation of the effective resistance on 
distance-regular networks [1], where the calculation is based on stratification introduced in [6] 
and Stieltjes transform of the spectral distribution (Stieltjes function) associated with the net- 
work. Also, in Ref. [1] it has been shown that the resistances between a node a and all nodes 
/3 belonging to the same stratum with respect to the a (-R^^ci) , /9 belonging to the i-th stratum 
with respect to the a) are the same and the analytical formulas have been given for two-point 
resistances R^ij(i),i — 1, 2, 3 in terms of the size of the network and corresponding intersection 
array without any need to know the spectrum of the pseudo inverse L~^. In the next work 
[2], the authors have used the algebraic structure of distance-regular networks (Bose-Mesner 
algebra) such as stratification and spectral techniques specially the well known Christoffel- 
Darboux identity [8] from the theory of orthogonal polynomials to give a recursive formula for 
calculation of all of the resistance distances Ra/sii), i — 1,2, ...,d {d is diameter of the network 
which is the same as the number of strata) on the network without any need to calculating the 
spectrum of the pseudo inverse of the Laplacian of the network denoted by L~^. In this way 
they have shown that, the effective resistance strictly increases by increasing the shortest path 
distance defined on the network, i.e., i?Q,^(m+i) — Rafsi^) > for all m = 1, 2, d — 1. Here in 
this work, evaluation of effective resistances on more general networks called pseudo-distance- 
regular networks [3] or QD type networks [4] is investigated, where we use the stratification of 
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these networks and show that the effective resistances between a given node such as a and all of 
the nodes /3 belonging to the same stratum with respect to a are the same. Then, based on the 
spectral techniques, an analytical formula for effective resistances Rai3{^) such that — L^^ 
(those nodes a, of the network such that the network is symmetric with respect to them) 
is given in terms of the first and second orthogonal polynomials associated with the network. 
Particularly, due to the fact that, in distance-regular networks we have = L^^ for all nodes 
a, /9 of the network, the effective resistances Rapim-) for m — l,2,...,d are calculated directly, 
by using the given formula. 

The organization of the paper is as follows. In section 2, we give some preliminaries such 
as definitions related to graphs, their adjacency matrices, stratification of the graphs, pseudo- 
distance-regular graphs and spectral distribution associated with the graphs. In section 3, an 
analytical formula for calculating the effective resistances in pseudo-distance-regular graphs as 
resistor networks is given in terms of the orthogonal polynomials of the first kind and second 
kind associated with the networks. The paper is ended with a brief conclusion and an appendix. 

2 Preliminaries 

In this section we give some preliminaries such as definitions related to graphs, corresponding 
stratification, pseudo-distance-regular graphs and spectral distribution techniques. 

2.1 Graph and its adjacency matrix 

In this section, we review the stratification of the graphs and the notion of pseudo-distance- 
regularity. 

A graph is a pair F = {V,E), where F is a non-empty set and £■ is a subset of {{x,y) : 
x,y E V,x y}. Elements of V and of E are called a vertex and an edge, respectively. Two 
vertices x,y E V are called adjacent if {x, y) e E, and in that case we write x ~ y. For a 
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graph r = (V, E), the adjacency matrix A is defined as 

r 1 if q; ~ /? 

{A)a,p = . (2-1) 

otherwise 

Conversely, for a non-empty set y, a graph structure is uniquely determined by such a matrix 
indexed by V. 

The degree or valency of a vertex x & V is defined by 

K(x)^\{yeV:yr^x}\ (2-2) 

where, | • | denotes the cardinahty. The graph is called regular if the degree of all of the vertices 
be the same. In this paper, we will assume that graphs under discussion are regular. A finite 
sequence xo,Xi, ...,Xn G V is called a walk of length n (or of n steps) if Xi^i ~ Xi for all 
i = 1,2, Let f{V) denote the Hilbert space of C-valued square-summable functions on 

V. With each (3 & V we associate a vector such that the /3-th entry of it is 1 and all of 
the other entries of it are zero. Then : e V} becomes a complete orthonormal basis of 
l'^{V). The adjacency matrix is considered as an operator acting in P{V) in such a way that 

m = E !«)• (2-3) 

2.2 stratification 

For X y let d{x, y) be the length of the shortest walk connecting x and y. By definition 
d{x, x) = for all x & V . The graph becomes a metric space with the distance function d. 
Note that d{x, y) = 1 if and only if x ~ y. We fix a point o e F as an origin of the graph, 
called reference vertex. Then, the graph F is stratified into a disjoint union of strata: 

oo 

y = U F,(o), F,(o) := {a e y : d{a, o) = (2-4) 

i=0 

Note that Fj(o) = may occur for some i > 1. In that case we have Fj(o) = Fj_|_i(o) = ... = 0. 
With each stratum Fi(o) we associate a unit vector in l'^(y) defined by 

\ E (2-5) 
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where, Ki = |ri(o)| is called the i-th valency of the graph := |{7 : 9(0,7) = 01 — |ri(o)|). 
The closed subspace of l'^{V) spanned by is denoted by A(r). Since becomes a 

complete orthonormal basis of A(r), we often write 



A(r) = $:©c|(/.,). 



(2-6) 



2.3 Pseudo-distance-regular graphs 

Given a vertex o; e ^ of a graph F, consider the stratification (2-4) with respect to a such that 
ri{a) — ior i > d. Then, we say that F is pseudo-distance- regular around vertex a whenever 
for any (3 e Tk{a) and < k < d, the numbers 



1 11 

^^^^^ ^ wv^ ^ ^^^^^ ^^TT^T ^ 

7eri(/3)nrfc_i(a) '^^^^^ 7eri(/3)nrfc(a) 7eri(/3)nr,+i(a) 

(2-7) 



do not depend on the considered vertex (3 e Ffe(Q;), but only on the value of k. In such a case, 
we denote them by Ck, Qfe and bk respectively. Then, the matrix 



V 



(2- 



ci ... Cd-i Cd 
flo ai ... ttd-i ttd 
bo bi ... bd-i 

is called the (pseudo-)intersection array around vertex a of F. It is shown in Ref. [3] that this 
is a generalization of the concept of distance-regularity around a vertex (which in turn is a 
generahzation of distance-regularity). It should be noticed that for regular graphs F {k,{/3) — k 
for all /3 e y), the numbers c^, and bk read as 



Cfc= |Fi(/3)nFfe_i(a)|, afe = |Fi(/?)nFfe(a)|, 6^ = |Fi(/?) n Ffc+i(a)|, 



(2-9) 



where we tacitly understand that F_i(q;) = (see Fig. 1). The intersection numbers (2-9) and 

the valencies = |Fi(a)| satisfy the following obvious conditions 



+ + Cj — AC, Ki-ibi-i — KiCi, i — 1, ...,d, 
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Ko = ci = 1, bo = Ki = K, (co = bd = 0). (2-10) 

One should notice that, the definition of pseudo-distance regular graphs together with the 
Eq.(2-10), imply that in general, the valencies ( the size of the i-th stratum) for i = 0, 1, o? 
do not depend on the considered vertex /3 e ^k{o(), but only on the value of k. 

The notion of pseudo-distance regularity has a close relation with the concept of QD type 
graphs introduced by Obata [4], such that for the adjacency matrices of this type of graphs, 
one can obtain a quantum decomposition associated with the stratification (2-4) as 

A^A+ + A- + A°, (2-11) 

where, the matrices A^, A~ and ^4° are defined as follows: for /3 e ^k{c(), we set 

Af)s if 5 e Tk+i{a) 



(^+) 



135 



otherwise 

Af3s if 5 e Tk_i{a) 

otherwise 
Aps if 5 e Tkia) 

otherwise 



Or equivalently, for 



A-^%(5)^ Y: \k+l,S), A-\k,P)^ J2 \k-l,S),, A'\k,P)^ Yl 

Serk+i{a),5r^P SeTk-i{a),Sr^f3 5erfc(a),5~/3 

(2-12) 

Since P e Fkia) and /3 ~ 5 then S e Tk-iia) U Fkia) U Tk+iia). 

It has been shown in Ref. [4] that, if A(r) is invariant under the quantum components 
A'^,A~ and A^, then there exist two sequences (called Szego- Jacobi sequences) {u!i}iZi and 
{aiyfZi derived from A such that 

A-|0o) = O, A-\<f>l)^^l\<^>l_^), 1>1, 
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A^\<Pi) = ai\<Pi),l>0, (2-13) 

where ui+i = ^/t^ (j), = \{i E Ti{a) : i ~ j}\ for j e Ti+i{a) and a; = ko(j), such 

that ACo(j) — |{^ G ^ ill for j e r/(Q;), for / > 0. One can easily check that the coefficients 
ctj and cui are given by 

ak = ak^ K-bk- Ck, LOk = 01^ h-iCk, k ^ 1, d. (2-14) 
By using (2-11) and (2-13), one can obtain 

A|0,)=A+i|0i+i)+«d0/)+A|0i-i), ^>0, (2-15) 
with Pi :— y/uTi. Then, by using (2-15), one can deduce that 

\<l>i)^Pi{A)\<l>o), l^l,2,...,d, (2-16) 

where, Pi{A) = 00 + 01^4 -|- 02^!^ + ... + aiA'' is a polynomial of degree / in indeterminate A and 
conversely ^4' can be written as a linear combination of Po{A), Pi{A), Pd{A). In fact, it can 
be shown that [9], the unit vectors \(f)i) for 1 — 0,1, d are the orthonormal basis produced 
by application of the orthonormalization process (Lanczos algorithm) to the so called Krylov 
basis {|0o), ^|0o), ^^|0o), A'^lcpo)} (for more details see for example [6], [9]). 

2.4 Spectral distribution associated with the network 

In this subsection we recall some facts about the spectral distribution associated with the 
adjacency matrix of the network. In fact, the spectral analysis of operators is an important 
issue in quantum mechanics, operator theory and mathematical physics [10, 11]. Since the 
advent of random matrix theory (RMT), there has been considerable interest in the statistical 
analysis of spectra [12, 13, 14]. Also, the two-point resistance has a probabilistic interpretation 
based on classical random walker walking on the network. Indeed, the connection between 
random walks and electric networks has been recognized for some time (sec e.g. [15, 16, 17] 
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), where one can establish a connection between the electrical concepts of current and voltage 
and corresponding descriptive quantities of random walks regarded as finite state Markov 
chains (for more details see [18]). Also, by adapting the random- walk dynamics and mean- 
field theory it has been studied that [19], how the growth of a conducting network, such as 
electrical or electronic circuits, interferes with the current flow through the underlying evolving 
graphs. In Ref. [1], it has been shown that, there is also connection between the mathematical 
techniques introduced in previous subsections and this subsection such as Hilbert space of 
the stratification and spectral techniques (which have been employed in [6, 7, 9, 20, 21] for 
investigating continuous time quantum walk on graphs), and electrical concept of resistance 
between two arbitrary nodes of regular networks, and so the same techniques can be used for 
calculating the resistance. Note that, although we take the spectral approach to define the 
effective resistance in terms of orthogonal polynomials (which are orthogonal with respect to 
the spectral distribution /i associated with the network) with three term recursion relations, 
in practice as it will be seen in section 3, the effective resistances will be calculated without 
any need to evaluate the spectral distribution /i. 

It is well known that, for any pair {A, |0o)) of a matrix A and a vector \(f)o), it can be 
assigned a measure /x as follows 

//(x) = (0o|£;(x)|</)o), (2-17) 

where E{x) = J2i \ui){ui\ is the operator of projection onto the eigenspace of A corresponding 
to eigenvalue x, i.e., 

A^ j xE{x)dx. (2-18) 
It is easy to see that, for any polynomial P{A) we have 

P{A) = j P(x)E(x)dx, (2-19) 

where for discrete spectrum the above integrals are replaced by summation. Therefore, using 
the relations (2-17) and (2-19), the expectation value of powers of adjacency matrix A over 
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starting site \4>q) can be written as 

((/.olA^lc/.o) = / xXda;), m = 0,1,2,.... (2-20) 

JR 

The existence of a spectral distribution satisfying (2-20) is a consequence of Hamburgers the- 
orem, see e.g., Shohat and Tamarkin [[22], Theorem 1.2]. 

Obviously relation (2-20) implies an isomorphism from the Hilbert space of the stratification 
onto the closed linear span of the orthogonal polynomials with respect to the measure /j,. More 
clearly, the orthonormality of the unit vectors implies that 

^ij = {^M = I Pi{x)Pj{x)ii{dx), (2-21) 

J R 

where, we have used the Eq. (2-16) to write \4>i) = Pi{A)\4>Q). Now, by substituting (2-16) in 
(2-15), we get three term recursion relations between polynomials Pj{A), which leads to the 
following three term recursion relations between polynomials Pj{x) 

xPk(x) = (3k+iPk+i{x) + akPk(x) + pkPk-i{x) (2-22) 

for k — 0, d — 1, with Po{x) — 1. Multiplying two sides of (2-22) by Pi...Pk we obtain 

(3^...(3kxPk{x) = Pi...Pk+iPk+i{x) + akPi...PkPk{x) + Pl.Pi...Pk-iPk-i{x). (2-23) 

By rescaling as Qk — Pi...PkPk, the spectral distribution /i under question is characterized 
by the property of orthonormal polynomials {Qk} defined recurrently by 

Qo{x)^l, Qi{x)^x, 

xQk (x) = Qk+1 (x) + akQk (x) + PlQk-i {x), k>l (2-24) 

(for more details see Refs.[8, 22, 23, 24]). 

It should be noticed that, the starting values of the recurrence (2-24) are Q_i = 0, Qo = 1- 

If one starts from q_i = —1, Qo = and uses the same recurrence (with uq = 1) 

Xqnix) = Qn+lix) + an9n(a;) + UnQn-lix), (2-25) 
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then Qn is of degree n—1, and by Favard's theorem the different g„'s are orthogonal with respect 
to some measure. The g„'s are called orthogonal polynomials of the second kind (sometimes 
for Qn we say that they are of the first kind) . They can also be written in the form 



(for more details see for example [25]). 

3 Two-point resistances in regular resistor networks 

A classic problem in electric circuit theory studied by numerous authors over many years, is 
the computation of the resistance between two nodes in a resistor network (see, e.g., [26]). 
The results obtained in this section show that, there is a close connection between the tech- 
niques introduced in section 2 such as Hilbert space of the stratification and the orthogonal 
polynomials of the first and second kind associated with the networks and electrical concept 
of resistance between two arbitrary nodes of the networks. 

For a given regular network F with v vertices and adjacency matrix A, let r^j — rji be the 
resistance of the resistor connecting vertices i and j. Hence, the conductance is Cij — r^^ — Cji 
so that Cij = if there is no resistor connecting i and j. Denote the electric potential at the 
i-th vertex by Vi and the net current flowing into the network at the i-th vertex by /j (which 
is zero if the i-th vertex is not connected to the external world). Since there exist no sinks 
or sources of current including the external world, we have the constraint J2^=i h = 0. The 
Kirchhoff law states 




(2-26) 



V 



i = 1,2,. ..,v. 



(3-27) 



Exphcitly, Eq.(3-27) reads 
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where, V and / are vectors whose components are Vi and li, respectively and 

L = E^^\^){^\-E^^J\^){J\ (3-29) 



is the Laplacian of the graph F with 



Ci= Cij, (3-30) 



for each vertex a. Hereafter, we will assume that all nonzero resistances are equal to 1, then 
the off-diagonal elements of —L are precisely those of A, i.e., 

L = kI -A, (3-31) 

with K = K,i = deg{a), for each vertex a. It should be noticed that, L has eigenvector 
(1, 1, 1)* with eigenvalue 0. Therefore, L is not invertible and so we define the psudo-inverse 
of L as 

= Y: K'E., (3-32) 

where, Ei is the operator of projection onto the eigenspace of L"^ corresponding to eigenvalue 
Aj. It has been shown that, the two-point resistances R^jj are given by 

R^p = {a\L-'^\a) + {P\L-^\P) - - (3-33) 

This formula may be formally derived using KirchofF 's laws, and seems to have been long 
known in the electrical engineering literature, with it appearing in several texts, such as Ref. 
[27]. 

Now, consider two nodes a,P such that L^^ = L^^ (as we will see in subsection 4.1, 
for distance-regular graphs as resistor networks, the diagonal entries of are independent 
of the vertex, i.e., for all a,P eV, we have — L'^g). Then, from the relation (3-33) and 
the fact that is a real matrix, we can obtain the two-point resistance between the nodes 
a and /3 as follows 

Rap = 2{L-'^-L-'g)- (3-34) 
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It should be noticed that, due to the stratification of the network, all of the nodes belonging 
to the same stratum with respect to the reference node (a node which stratification is done 
with respect to it), possess the same effective resistance with respect to the reference node 
(the proof is given in the appendix) . More clearly, in order to evaluate the effective resistance 
between two nodes a and /3 of a network, we consider one of the nodes, say a, as reference 
node and stratify the network with respect to a. Then, j3 will belong to one of the strata with 
respect to a, say the m-th stratum r^(Q;). Now, for all /? e r^{a), one can write 



--1 



= {a\L-'\(3) = -L(a|L-Vm) = -L(a|P^(^)L-». (3-35) 

a/ Km. \/ Km 



Jr-\k\ K — X 



where, we have used the equations (2-5) and (2-16) and the lemma given in the appendix. 
Then, by using (3-34), we obtain for all (3 e Tm{oi) 

D ^ r / r-l I TD ( A\T-l\ 1 ^ / I V l^m^l---^mX — Qm\A) , , 

Rapim) = ^={V^^aa - {Pm{A)L )„„) = , {a\ ^ \a} = 

Urn s/ l^m^l---^m ^ 

\J l^m^l---^m JR-{k} 

where, the upper index m in L~^(m) and Ra/3(m) indicate that (3 belongs to the m-th stratum 
with respect to a. Note that, we have substituted Pm{x) — , Qm{x^ and used the equality 
Qm{i^) = i/Km'-^i---'^m which Can be verified easily. It should be noticed that, the result (3-36) 
can be written as 

Now, by using (2-26) and (3-37), we obtain the main result of the paper as follows 

2 Id 

Raf3(m) = , :{gm(^) " ^{^Qm{x))\^=^} , m = 1, 2, d. (3-38) 

We recall that, the result (3-38) can be used only for evaluating the effective resistance between 
nodes such that the network is symmetric with respect to them in the sense that, stratification 
of the network with respect to these nodes produces the same strata. Then for such nodes a 
and P, we will have = L^^. 
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4 Examples 

4.1 Evaluation of effective resistances on distance- regular networks 

First, we recall the definition of special kind of pseudo-distance-regular networks called distance- 
regular networks: 

Definition. A pseudo-distance regular network F = (K E) is called distance-regular with 
diameter d if for all k e {0, 1, ■■■,d}, and a,P & V with (3 e Ffe(Q;), the numbers Ck{(3), afe(/9) 
and bk{/3) defined in (2-7) depend only on k but do not depend on the choice of a and /3. 

It should be noticed that, in distance-regular networks, the i-th adjacency matrix of the 
network T — {V,R) is defined by 

f 1 if d(a,(3) =i 
(A)a,/3 = . (4-39) 

I otherwise 

Then, from the definition of A^, for the reference state |0o) (|0o) — \o), with o eV as reference 
vertex), we have 

A|0o)= E 1/5)- (4-40) 
Then by using (2-5) and (4-40), we have 

A|0o) = y/^M- (4-41) 

Also, it can be shown that, for adjacency matrices of a distance regular graph, we have 

AiAi = bi_iAi_i + ttiAi + Q+iA+i, for i = l,2,...,d- 1, 

AiAd = bd-iAd-i + aaAd. (4-42) 

Using the recursion relations (4-42), one can show that A^ is a polynomial in Ai of degree i, 
i.e., we have 

Ai^PiiAi), i = l,2,...,d, (4-43) 
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and conversely A\ can be written as a linear combination of /, ^41, vl^ (for more details see 
for example [6]). 

Now, it should be noticed that, stratification of distance- regular networks will be indepen- 
dent of the choice of the reference node (the node which stratification is done with respect to 
it). Then, clearly we will have — L^^ for all a, P E V with a P and consequently, the 
result (3-38) can be used for evaluation of the effective resistance between any two arbitrary 
nodes a, It is sufficient to choose one of these nodes, say a, as reference node and stratify the 
network with respect to it. Then, /3 will be contained in one of the strata, say m-th stratum, 
with respect to a and so the effective resistance between a,(3 E r^(Q;) can be evaluated via 
Eq.(3-38). 

1. Cycle network Cy 

A well known example of distance-regular networks, is the cycle network Cy with k = 2. The 
network C„ ior v — 2k or v — 2k + 1 consists of A; -|- 1 strata, where the intersection arrays for 
even and odd number of vertices are given by 

{bo, ci, = {2, 1, 1, 1; 1, 1, 2} (4-44) 

and {6o, ...,6fe-i;ci,...,Cfc} = {2, 1,...,1;1,...,1,1}, (4-45) 
respectively. Then, by using (2-14), for even v — 2k the QD parameters are given by 

ai — 0, i — 0,1, ...,k; uJi—uJh — 2, uji — l, i — 2,...,k — l, (4-46) 

where, for odd v — 2k + 1, we obtain 

ai — 0, i — 0,1, ...,k — 1, ak — 1; uJi — 2, uji — l, i — 2,...,k. (4-47) 

We consider v — 2k (the case v — 2k + 1 can be considered similarly). Then, by using (2-24) 
and (4-46), one can obtain 

Qo = l, Qi = x, Q,{x) = 2Ti{x/2), i > 1; q,{x) = Q,_^{x), i>l, (4-48) 
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where, Tj's are Chebyshev polynomials in one variable, which are recursively defined by 

To = l, Ti=x, T„(a;) = 2a;T„_i - T„_2, n > 1. (4-49) 

Then, by using (3-38) and (4-48), the effective resistance between any two nodes a, (3 E Tm{Q:), 
is obtained as 

Id 1 
Rapm = gi(2) - ^-Q^Qii^)\-=2 = 1 - ^ ' 

= 2r^_i(l) - ^r;(l), for m > 1. (4-50) 

2. (d-|-l)-Simplex fractals 

(d-l-l)-simplex fractal is a generalization of a two dimensional Sierpinski gasket to d-dimensions 
such that its subfractals are {d + l)-simplices or d-dimensional polyhedra with S^^d + 1)- 
symmetry. In order to obtain a fractal with decimation number 2, we choose a (ci-|- 1) -simplex 
and divide all the links (that is the lines connecting sites ) into 2 parts and then draw all pos- 
sible d-dimensional hypcrplanes through the links parallel to the transverse d-simplices. Next, 
having omitted every other innerpolyhedra, we repeat this process for the remaining simplices 
or for the subfractals of next higher generation. This way through (d + l)-simplex fractals are 
constructed. We label subfractals of generation (n -|- 1) in terms of partition of 1 into {d + 1) 
positive integers Ai, A2, A^+i. Each partition represents a subfractal of generation n, and A 
shows the distance of the corresponding subfractal from d-dimensional hyperplanes which con- 
struct the ((i-l-l)-simplex. On the other hand, each vertex denoted by partition of 2 into (d+l) 
non-negative integers 771, 772, , Vd+i ^-nd obviously the i-th vertex of subfractal (Ai, A2, ^d+i) 
is denoted by rjj = Xj + Sij, where j = 1, 2, , 0? -|- 1 (for more details see [28], [29]). 

Now, we consider the {d + l)-simplex fractal with decimation number 6 = 2 such that all 
of the d+1 vertices (20. ..0), (020. ..0), (0...02) are connected via a resistor to each other 

d+l d+1 d+1 

(see Fig. 2 (a) and Fig. 2 (b) for 0? = 2 and d = 4, respectively). Then, the number of vertices 
of {d + l)-simplex fractal is = Cf+^ + C^^ = {d+W+2) ^^^^ ^^^^ ^^^^^^ ^^-^ ^^^^^^ 

is K, — 2d. Also, it can be easily shown that the network has 3 strata with respect to the 
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reference node (200. ..0), where the unit vectors i = 0, 1, 2 are given as follows 

|0o) = |2OO^, 

d+l 

= ^^(|00...0 2 0...0) + I10...0 1 0...0)), 



l« = \/s(d^_S |0...0^0...0^0...0). (4-51) 
Then, one can show that 

A|0o) = V2d|0i), 



A|02) = 2v/d^|0i)+2(d -2)102). (4-52) 

Also, by using the recursion relations (4-52) one can easily show that the adjacency matrices 
A = Ai, A2 — J — A — I (where, J is all one matrix) satisfy the following relations 

A' = 2d.I (d+l)(d+2) + d.A + 4^2, 

2 

AA2^ {d-l)A + 2{d-2)A2. (4-53) 

It should be noticed that if we stratify the network with respect to another reference node 
such as (1100. ..0), the unit vectors will be obtained as 




1) = (I10...0 1 0...0) + 1010...0 1 0...0)) + 1200...0) + 1020...0)}, 

y2d j=3,...,d+i ^ ^ 



102)- W J E |00...0^0...0)+ J2 |00...04_^0...04_^0...0)}. (4-54) 

V V / i=3,...,d+l j i,j=3,...,d+l;i<j j j 

Then, one can show that, the same recursion relations as in (4-52) and (4-53) are satisfied 
for this stratification, i.e., stratification of the network gives three-term recursion relations 
independent of the choice of the reference node and so the network is distance-regular. 
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Now, by using the equations (2-15) and (4-52), one can obtain 

ao = 0, ai = d, a2 = 2{d - 2) ; uui = 2d, uj2 = 4((i - 1). (4-55) 

Then, by using the recursion relations (2-24) and (2-25), one can see that 

Qi{x) = X, Q2{x) = — dx — 2d ; qi{x) = 1, q2{x) = x — d. (4-56) 

It should be noticed that, the intersection array of the {d+ 1)- simplex fractal with decimation 
number b — 2, can be evaluated by using (2-14) and (4-55) as: {bo, bi, ci, 02} — {2d, d — l;l, 4} 
which imphes that k ^ bo ^ 2d, K2 = ^ = Now, by using (3-38) and (4-56), the 

effective resistance between any node a & V and /3 e r^{a) for m = 1, 2 is given by 

R = 2 J_ 1 _ 2 ^ ^ + 3 

^a,pm ^{ ^) ^\ (d+i)(d + 2)^ (d+l)(d + 2)' 

^-^^^^ ^ d(d^^^^('^)-^('"-^)l^-'^> = d{hl)^^- {d+l){d + 2f ^^ = {dA)il+2) 

(4-57) 

3. Hexagon network 

Consider the hexagon network with its diameters shown in Fig. 3. As the figure shows, this 
network has 6 nodes with intersection array {60, ^1; Ci, C2} = {3, 2; 1, 3}. Then by using (2-10) 
and (2-14), one can obtain 

K — 3, K2 — 2; ao — ai — a2 — 0, cui — 3, L02 — 6. (4-58) 

Then, by using (2-24) and (4-58), we obtain 

<3o = 1, Q\^x, Q2{x)^x^-'i ; gi(x) = 1, 52(2;) ^x. (4-59) 

Now, by using (3-38) and (4-59), the effective resistance between any node a & V and /? e 
r^{a) for m = 1, 2 is given by 

Ra,pW = ^(1 -l) = l^ ^a,/J(^) = IMS) - 1{2X)U=3} = | (4-60) 
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4. A bipartite distance-regular network with 2n nodes 

Consider a distance-regular network with 2n nodes and adjacency matrices as 

= (7^0 (J„ -4), h^iJn- In), A3 = (7^ ® 7„, (4-61) 

where, J„ is an n x n matrix such that all of its entries are one. Then, one can show that 

Al^{n-l)l2n + {n-2)A2, A1A2 ^ {n - 2)Ai + {n - l)As, AiAs ^ A2, (4-62) 

By using (4-42) and (4-62), the intersection array of the network is given by 

{60, bi, 62; ci, C2, C3} = {n - 1, n - 2, 1; 1, n - 2, n - 1}. 

Then, by using (2-10) and (2-14), one can obtain 

K — n — 1, K2 — n — 1, Ks — 1, 

ai — 0, i = 0, 1, 2, ; ui — n — 1, uj2 — {n — 2)^, uj^ — n — 1. (4-63) 
Now, by using (2-24) and (4-63), we obtain 

qi{x)^l, q2{x)^x, q^{x)^x'' -{n-2f. (4-64) 

So, by using (3-38) and (4-64), the effective resistance between any node a & V and /3 e Tmicy) 
for m = 1, 2, 3 is given by 

2 1 2n - 1 

2 , , , 1 , , 2(n - 1) 

^-^'^^ = (n-l)(n-2) ^^^^" - - ^(2.:)U^.-i} = 

2 1 2n — 3 

^-^'^^ = {n -l){n- 2) ^'^^^^ 2;^^^^' - ~ + ^> = (n-l)(n-2) - ^^"^^^ 
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4.2 Evaluation of effective resistance in examples of pseudo-distance- 
regular networks 

1. Pseudo-distance- reguleir network derived from Hadanicird network with 16 

nodes 

Consider the pseudo-distance-regular network shown in Fig. 4. This network is obtained from 
the Hadamard network with intersection array {4, 3, 2, 1; 1, 2, 3, 4}. As Fig. 4 shows, the net- 
work is symmetric with respect to the initial and final (horizontal) nodes 1 and 16 e r4(l) 
and also with respect to the initial and final (vertical) nodes 6 and 11 e r4(6). One should 
notice that stratification of the network with respect to the nodes 1 and 16 produces the same 
strata. For stratification with respect to the node 1 or 16, we have 

K = 4, K,2 = Q, 1^3 = 4, K.4 = 1, 

ai — 0, i — 0,...,3; uji — 4, UJ2 — UJ3 — 6, 004 — A. 
(see Eq. (2-13)). Then, by using the recursion relations (2-24) and (2-25), one can obtain 

Q^ix) ^ -16x^ + 24, q^ix) ^ x^ -12x. (4-66) 

Now, from Eq. (3-38), the effective resistance between nodes 1 and 16 e r4(l) is given by 

= ^{94(4) - ^(4x3 - 32x)U=4} = l- (4-67) 

In order to evaluate the effective resistance between vertical nodes 6 and 11, one must consider 
the stratification of the network with respect to the node 6 or 11. For this stratification, we 
obtain 

K — A, K2 — 4:, K3 — 4, K4 — 3. 
ai — 0, i — 0,...,3; cui — A, 002 — A, cus — l, 004 — 12. 
Then, by using the recursion relations (2-24) and (2-25), one can obtain 

Q4ix) = x"^ - 9x^ + A, q4{x)^x^-5x. (4-68) 
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Therefore, from Eq. (3-38), the effective resistance between nodes 6 and 11 e r4(6) is given 

by 

Re,.. = Y^{94(4) - -{4x' - 18x)U=4} = ^. (4-69) 

2. Pseudo-distance-reguleir network derived from DesEirgues 

This network has 20 nodes with 

{bo, bi, 62, &3, &4; ci, C2, C3, C4, C5} = {3, 2, 2, 1, 1; 1, 1, 2, 2, 3}. 
Then, by using (2-10) and (2-14), one can obtain 

K = 3, K2 = 6, = 6, ft;4 = 3, ^5 = 1 

ctj = 0, i = 0, 1,...,5; LUi = 3, a;2 = 2, a;3 = 4, a;4 = 2, 00^ = 3. 

The stratification with respect to the initial node 1 and the final node 20 produces the same 
strata. Therefore, the effective resistance between the nodes 1 and 20 can be evaluated by 
using the Eq. (3-38). By using the recursion relations (2-24) and (2-25), we obtain 

Q^ix) ^ -nx^ + 22x, 55(2;) = - 8x^ + 4. (4-70) 

Then, from Eq. (3-38), the effective resistance between nodes 1 and 20 e r5(l) is given by 

^..0 = ^{?5(3) - ^{5x' - 33x' + 22)U=3} = (4-71) 

3. Pseudo-distance-reguleir network derived from Hadameird network with 32 

nodes 

Consider the pseudo-distance-regular network obtained from the Hadamard network with in- 
tersection array {8,7,4, 1; 1,4,7,8} such that, the network is symmetric with respect to the 
initial and final (horizontal) nodes 1 and 32 e r4(l) and also with respect to the initial and 
final (vertical) nodes 10 and 23 e r4(10). One should notice that stratification of the network 
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with respect to the nodes 1 and 32 produces the same strata. For stratification with respect 
to the node 1 or 32, we have 

K — 8, K2 — 14, ^3 = 8, K4 — 1, 

ttj = 0, i = 0, 4; LUi = 8, UJ2 = 28, cu^ = 28, cu^ = 8. 
(see Eq. (2-13)). Then, by using the recursion relations (2-24) and (2-25), one can obtain 

Q4(^) =2:^-64x2 + 224, q^{x) ^ - 56x. (4-72) 

Now, from Eq. (3-38), the effective resistance between nodes 1 and 32 e r4(l) is given by 

^1.32 = ^2^^'^^^ ~ h^^""' ~ 128a;)U=8} = f (4-73) 

In order to evaluate the effective resistance between vertical nodes 10 and 23 e r4(10), one 
must consider the stratification of the network with respect to the node 10 or 23 and evaluate 
R-^Q 23 before. 
5. A network with 9 nodes 

Consider the network given in Fig. 5 with 9 nodes and the following intersection array 

{&o,&i;ci,C2} = {4,2;l,2}. 

Then by using (2-10) and (2-14), one can obtain 

K = 4, = 4, ; cto = 0, «! = 1, a2 = 2, a;i = a;2 = 4. 

As it can be seen from Fig. 7, the stratification with respect to the nodes 1 and 9 e r2(l) 
produces the same strata. Then, by using the recursion relations (2-24) and (2-25), one can 
obtain 

Q2{x) = x"^ - X - 4, q2{x) = X -1. (4-74) 
Now, from Eq. (3-38), the effective resistance between nodes 1 and 9 e r2(l) is given by 

^,9 = ^{?2(4) - l{2x - 1)U=4} = l- (4-75) 
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6. A network with 16 nodes 

Consider the network given in Fig. 6 with 16 nodes and the following intersection array 

{bo, bi, 62, fes; ci, 02, C3, C4} = {4, 3, 2, 1; 1, 2, 3, 4}. 

Then by using (2-10) and (2-14), one can obtain 

K — A, K2 — 6, K3 — 4, K4 — 1; ai — 0, i — 0, 1, 4, uji — 4, 002 — — 6, 004 — 4. 

As it can be seen from Fig. 8, the stratification with respect to the nodes 1 and 16 e r2(l) 
produces the same strata. Then, by using the recursion relations (2-24) and (2-25), one can 
obtain 

Q2ix) q2{x) = X. (4-76) 

Now, from Eq. (3-38), the effective resistance between nodes 1 and 16 e r2(l) is given by 

^M6 = - ^(2^)U=4} = ^. (4-77) 

7. Generalized G2{m) type network 

Consider the network with 2(2'""'"^ — 1) nodes and the following intersection array 

{60, bi, bm] Cl, C2, ...Cm, Cm+l} = {3, 2, 2; 1, 1, ...1, 3}. 
Then by using (2-10) and (2-14), one can obtain 

K, = 3, Hi = 3.2 , i = 2, m, ACm+i ~ 2"^, 

ctj = 0, i = 0,1, ...,m + 1; a;i = 3, U2 = ... = uJm = ^m+i = 6. 

Then, the stratification with respect to the initial node 1 and the final node 2"*+^ — 2 e ri(l) 

produces the same strata. Therefore, from the fact that Qi{x) = x and qi{x) = 1, the effective 

resistance between nodes 1 and 2"*+^ — 2 e ri(l) is obtained as 

2 1 2"*+^ — 3 

-^i,2m+2_2 ~ 3^"*" ~ 2(2"*+^ — 1)'^ ~ 3('2m+i _ ■ (4-78) 
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Based on the stratification of tlie pseudo-distance-regular networks and using spectral tech- 
niques, the evaluation of the effective resistances on these networks was discussed. It was shown 
that, in these types of networks, the effective resistances between a node a and all nodes /3 
belonging to the same stratum with respect to the a are the same. Then, an explicit analyt- 
ical formula for the effective resistance between two nodes a, /3 of a pseudo-distance-regular 
resistor network such that = L^^ iL~^ is the pseudo-inverse of the Laplacian matrix of 
the network) was given in terms of the first and second orthogonal polynomials associated 
with the network. It was deduced that, the obtained result can be used for evaluation of the 
effective resistance between any two arbitrary nodes a, (3 in distance-regular networks, where 
we have — Lgl for all nodes a, p. 
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Appendix 

In this appendix we prove the following lemma in connection with the equality of effective 
resistance between the reference node and all of the nodes belonging to the same stratum with 
respect to the reference node. 

Lemma. Let Rap denote the effective resistance between nodes a,P eV. Then for pseudo- 
distance-regular networks, by choosing one of the nodes, say a as reference node, the effective 
resistance Ra/s is the same for all nodes P e r^(a), where m e {1, 2, d}. 
proof. 

In order to prove the above lemma, we prove that in general for any pseudo- distance-regular 
network with diameter d, we have {(f)o\f{A)\P) (|0o) = \<^)) is the same for all P e ^i{oi), i.e., 

{^,\f{Am^^{<f>o\fiA)\<Pi), (A-i) 

where, f{A) is any function of the adjacency matrix A of the network such that f{A) — 
EtoOfc^*^ and 

^ Yl U)' l = OA,-,d. (A-ii) 



To this aim, we take the Fourier transform of the unit vectors for / = 0, 1, o? as follows 

— i: 



\<f>i,k)^^ E e'-'^'^/^'lj), l = 0,l,...,d; fc = 0,l,...,«;;-L (A-iii) 



Now, we show that 

((/)O|/(A)|0,,O)7^O, 

{Mf{AMi,k) = o, y k = i,...,Ki-i. (A-iv) 

To do so, we use the fact that there is a correspondence between the basis I,A,..., A'^~^ and the 
orthogonal polynomials Pi{A) defined by Eq.(4-43). In fact, as it was regarded previously (see 
arguments about the Eq.(4-43) ), A' for I — 0, 1, d can be written as a hnear combination 

oiPo{A),P,{A),...,Pa{A). 
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It should be noticed that, in Krylov subspace projection methods, approximations to the 
desired eigenpairs of an n x n matrix A are extracted from a d-dimensional Krylov subspace 

K,{\cPo),A) = span{\cl>o),A\cj>o),---,A''-'\M, (A-v) 

where, |0o) is often a randomly chosen starting vector called reference state and d <^ n, i.e., 
the vectors \4>q) , A\4>q) , A'^~^\4>q) constitute a basis for the Krylov subspace Kd{\4>Q),A). 
Then, the application of the orthonormalization process (the Lanczos algorithm which is a 
modified version of the classical Gram-Schmidt orthogonahzation process) to the Krylov ba- 
sis {^'^|0o)}fc=o is equivalent to the construction of a sequence of orthonormal basis — 
Pj{A)\(f)o), where Pj{A) — ao + aiA + ... + ajA^ is a polynomial of degree j in indeterminate A. 

As regards the above arguments, any function f{A) can be expanded as a linear combination 
of the polynomials Pj{A), i.e.. 



f{A)^^bjPj{A). 

j=0 



Then, we have 



d d 

/>o|/(A)|0;,fc) = J2bj{<Po\PjiA)\<pi,,) = J2b, {<Pj,o\<Pi,k) = 0, V A; = 1, - 1. 

j=0 



(A-vi) 



(A-vii) 



3=0 



Now, let we denote {(j)o\f{A)\j) by Xj for j G Ti{a). Then, from (A-vii) we have 



Xi 



X2 



^ (0o|/(^)|« ^ 




J \ y 

where, F is the ki x ki discrete Fourier transformation matrix(DFT) defined as Fj^ 
Therefore, by inverting F in (A-viii), we obtain 



(A-viii) 



1 ^2mjklKi 



Xi 



X2 



\ •'^1^1 J 



Ft 



^ (0o|/(^)|« ^ 




(0o|/(A)|0/,o) 



(A-ix) 
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That is, we obtain Xj = {(j)o\f{A)\j) = ;^(0o|/(^)|0i,o) for all j e Ti{a). 
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Figure Captions 

Figure-1: Shows edges through a and /3 in a pseudo-distance-regular graph. 
Figure-2: (a) Shows 3-simplex fractal with decimation number b — 2 (b) Shows 5-simplex 
fractal with decimation number 6 = 2. 

Figure-3: Shows the Hexagon network. 

Figure-4: Shows a pseudo-distance-regular network with 16 nodes derived from Hadamard 
graph. 

Figure-5: Shows a pseudo-distance-regular network with 9 nodes. 
Figure-6: Shows a pseudo-distance-regular network with 16 nodes. 
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